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ABSTRACT.  The  integral  of  the  elliptic  Gaussian 
HTitrTBution  over  a  centered  ellipse  may  be  evaluated 
after  simple  preliminary  transformations  from  tables 
of  the  circular  coverage  function.  The  probable  value 
of  the  distribution  over  a  centered  but  randomly  ori¬ 
ented  ellipse  where  all  orientations  are  equally  prob¬ 
able  may  be  evaluated  from  an  integral  involving  the 
circular  coverage  function. 
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FOREWORD 

The  estimation  of  the  probability  of  missiles  falling  within  speci¬ 
fied  areas  is  important  in  their  development  and  application.  For 
certain  applications  the  accurate  evaluation  of  this  probability  is 
required. 

This  report  presents  some  techniques  for  accurately  evaluating 
from  tables  of  the  circular  coverage  function,  the  probability  of  a 
missile  hit  within  a  centered  elliptical  area  under  various  assumed 
conditions  where  the  missile  deviations  are  normally  distributed. 

The  work  was  performed  under  Bureau  of  Naval  Weapons  Task 
Assignment  RUTO-3H  133/216  1/F008-03-001 ,  during  Fiscal  Year  1961. 

Review  of  this  report  for  technical  adequacy  was  performed  by 
D.  L.  Elliott  and  R.  E.  Wagner  of  this  Station. 


D.  J.  WILCOX,  Head 
Underwater  Ordnance 
^Department 


Released  under 
the  authority  of: 

WM.  B.  McLEAN 
Technical  Director 


NOTS  Technical  Publication  2766 
NAVWEPS  Report  7782 


Published  by . 

Manuscript . 

Collation . 

First  printing . 

Security  classification 


Underwater  Ordnance  Department 

.  807/MS- 138 

Cover,  10  leaves,  abstract  cards 

. 115  numbered  copies 

.  UNCLASSIFIED 


li 


1  Allegany  Ballistics  Laboratory,  Cumberland,  Md.  (Librarian) 

1  Applied  Physics  Laboratory,  University  of  Washington,  Seattle 
1  Convair  Scientific  Research  Laboratory,  San  Diego  (A.  L.  Berlad) 
1  Electric  Boat  Division,  General  Dynamics  Corporation,  Groton, 
Conn.  (Mr.  Atkinson) 

1  General  Electric  Company,  Defense  Electronics  Division, 
Pittsfield,  Mass.  (K.  I.  Igler) 

1  Hudson  Laboratories,  Columbia  University,  Dobbs  Ferry,  N.  Y. 

1  Ordnance  Research  Laboratory,  Pennsylvania  State  University 
(Development  Contract  Administrator) 

1  Scripps  Institution  of  Oceanography,  University  of  California, 

La  Jolla  (Document  Control) 

1  The  Rand  Corporation,  Santa  Monica,  Calif. 

1  Westinghouse  Research  Laboratories,  Pittsburgh  (Arthur  Nelkin) 
1  Woods  Hole  Oceanographic  Institution,  Woods  Hole,  Mass. 


NAVWEPS  REPORT  7782 


CONTENTS 

Introduction .  1 

Discussion  .  1 

Circular  Coverage  Function .  2 

Integral  of  the  Circular  Gaussian  Distribution  Over  a 

Centered  Ellipse  .  3 

Integral  of  the  Elliptic  Gaussian  Distribution  Over  a 

Centered  Circle .  7 

Integral  of  the  Elliptic  Gaussian  Distribution  Over  a 

Centered  Ellipse  .  9 

Integral  of  the  Elliptic  Gaussian  Distribution  Over  a 

Centered  Randomly  Oriented  Ellipse .  9 

Conclusions .  13 

Appendix:  Probability  Density  of  a  Randomly  Oriented  Elliptic 

Gaussian  Distribution . . .  14 

Nomenclature .  16 


References 


17 


NAVWEPS  REPORT  7782 


INTRODUCTION 

The  estimation  of  the  probability  of  a  missile  falling  within  a  speci¬ 
fied  area  under  various  assumed  conditions  has  been  the  subject  of 
study  for  many  years.  In  general,  it  is  assumed  that  the  range  and 
lateral  deviations  of  the  missile  are  normally  distributed  (i.  e.  , 

Gaus  sian-di  stributed) . 

Simple  cases  are  those  of  a  rectangular  target  area  with  sides  par¬ 
allel  to  range  and  lateral  dispersions  (assumed  uncorrelated),  and  a 
circular  target  area  centered  with  respect  to  a  circular  Gaussian  dis¬ 
tribution.  Recently,  other  cases  such  as  the  estimation  of  the  proba¬ 
bility  of  a  missile  falling  within  a  circular  target  area  where  the  dis¬ 
tribution  in  missile  deviations  is  offset,  circular  or  elliptical  Gaussian, 
has  received  study,  especially  by  the  Rand  Corp.  ,  Santa  Monica, 
California  (Ref.  1  to  12). 

In  general,  two  methods  are  available  for  estimating  the  probability 
of  a  missile  hit.  One  method  is  graphical  and  consists  in  plotting  the 
suitably  transformed  target  area  on  circular  probability  paper  and 
summing  the  probability  values  of  the  elements  within  the  target  area. 
The  other  method  consists  in  evaluating  the  probability  from  tabulated 
functions.  The  relative  convenience  of  the  methods  will  depend  upon 
the  nature  of  the  problem  involved.  However,  it  should  be  pointed  out 
that  summing  probability  elements  can  be  fatiguing  and  the  probability 
of  a  missile  hit  can  be  evaluated  more  accurately  from  tabulated  func¬ 
tions.  In  many  cases  such  accuracy  is  not  required,  but  in  cases 
where  the  number  of  missiles  thrown  is  large,  and  the  individual  proba¬ 
bilities  of  success  are  small,  the  accuracy  of  tabulated  functions  may 
be  required. 


DISCUSSION 

In  this  report  some  techniques  involving  the  circular  coverage 
function  are  discussed  for  estimating  the  probability  of  a  hit  in  the  use 
of  missiles  against  a  target  under  the  following  assumptions: 

1.  The  missile  range  and  lateral  deviations  are  normally  distributed 
and  known. 

2.  The  target  is  a  circle  or  an  ellipse. 

3.  The  center  of  the  target  and  center  of  impact  are  coincident. 

4.  The  orientation  of  the  target  with  respect  to  the  distribution  is 
known  or  it  is  random,  all  orientations  being  equally  probable. 
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In  order  to  avoid  unnecessary  discussion  of  well-known  familiar 
transformations,  reduction  to  standard  form  (i.  e.  ,  normalization)  will 
be  assumed  wherever  feasible. 


CIRCULAR  COVERAGE  FUNCTION 

The  circular  coverage  function  (Ref.  8)  is  the  probability  thatamissile 
will  hit  a  circle  of  radius  Rif  it  is  aimed  at  a  point  a  distance  r  from  the 
center  of  the  circle  and  if  it  is  subject  to  a  circular  Gaussian  impact- 
probability  law  of  unit  standard  deviation.  Reference  8  does  not  give  a 
derivation  of  this  function,  but  it  may  be  derived  as  follows. 

In  Fig.  1,  let  the  center  of  impact  be  represented  by  the  origin  O,  and 
let  the  target  be  represented  by  the  circle  with  center  C  and  radius  R. 
Without  loss  in  generality  an  x,  y- coordinate  system  with  origin  at  O  such 
that  C  lies  on  the  positive  x-axis  may  be  assumed.  Then  the  probability 
p(R,  r)  that  the  missile  will  hit  the  circle  is  given  by 


(1) 


p(R,  r)  =  — -  ff  e["(x2+y2)/2ldxdy 


2n  <x-r)r+?<R2 


Assuming  polar  coordinates  t,  0  with  the  pole  at  C  and  the  polar  axis 
along  the  positive  direction  of  the  x-axis 

(2)  x  *  r  +  t  cos  0 

y  =  t  sin  0 


and  Eq.  1  becomes 

(3)  p(R,  r)  *  —  eHr2/2)l  f  e[-(t2/2)lt  f  e("rt  cos  0)  d0  dt 

2  it  •'o  it 

From  Bessel  function  theory  (Ref.  13) 


(4)  e('rt  CO80)  *  IQ(rt)  +  2  £  (- l)nIn(rt)  cosn0 

n»l 

Substituting  Eq.  4  into  Eq.  3  and  integrating 

(5)  p(R,  r)  *  J'(r2/2)1  fR  e*'(t2/2)i0(rt)tdt 

~ o 

Equation  5  is  the  circular  coverage  function.  Certain  special  cases 
and  properties  of  this  function  should  be  noted. 

(6)  £(R.  0)  =  1  -  J'(r2/2)1 


Circular  Gaussian  probability  distribution 
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(7)  p(R,  «^  =  0  Target  infinitely  far  from  center  of  impact 

(8)  p(0,  r)  =  0  Target  of  null  area 

(9)  P(°°.  r)  =  1  Target  of  infinite  area 

(10)  p(-a,  P)  =  p(a,  -P)  =  p(-a ,  -p)  «  p(a,  P) 

Equations  6,  7,  8,  and  10  follow  directly  from  Eq.  5.  Equation  9  may 
be  obtained  by  substituting  (Ref.  13) 

(rt)2k 
22k(k!)2 

in  the  integrand  of  Eq.  5  and  integrating  term  by  term. 

Tables  of  the  circular  coverage  function  are  given  in  Ref.  8.  Ex¬ 
tensive  tables  (Ref .  9)  and  abridged  tables  (Ref.  11)  are  available  in 
which  the  function 


(12)  q(a,  p)  a  eHa2/2)]^  e[-(V2/2)]j  0(aV)VdV 

is  tabulated.  Since  from  Eq.  5,  9,  and  12 

(13)  p(R,  r)  =  1  -  q(R,  r) 

these  tables  may  be  used  to  obtain  p(R,  r).  The  function  q(R,  r)  is  the 
probability  that  the  missile  will  noFhit  the  target.  ~ 


INTEGRAL  OF  THE  CIRCULAR  GAUSSIAN 
DISTRIBUTION  OVER  A  CENTERED  ELLIPSE 


The  general  form  of  this  integral  is  given  by 


(14) 


p(a,  b)  *  — 


SL 


2 it  (x2/a2  +  y2/b2)  <  1 


t-(x2+y2)/2]  dxdy 
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In  what  follows  the  elliptic  semiaxes  a  and  b  are  assigned  their  positive 
values.  Under  the  transformation 

( 15)  x  *  t  cos  0 

y  *  t  sin  0 


Equation  14  becomes 

(16) 

where 
(17) 


p(a,  b)  a—  fZ  ft(9)  e[-<t2/2)]tdtde 
2 tr  Jo  Jo 


ab 


t(0)  = 


a2  s?n20  +  b2  cos20 


Integrating  Eq.  16  with  respect  to  t  and  substituting  Eq.  17  for  the  upper 
limit  of  integration 


(18) 


*  /*2w  g[-(  l/2)a2b2/(a2sin20+b2cos20)]  dg 

2tr  Jo 


2  r/Z  e[-(l/2)a2b2/(a2sin20+b2cos20)]  dQ 
n  Jo 


due  to  the  periodicity  of  the  integrand.  Applying  the  transformation 
(19) 

to  Eq.  18, 


b  t 

tan  0  -  —  tan  — 
a  2 


(20)  p(a.b).  1-  -  eMa2+bty4]  f"  e(-[(a2-b  2)/4]co.;| - 1 

ir  Jo  +  b2  + 


2abd& 


(a2  -  b^Jcos  £ 


In  order  to  evaluate  Eq.  20,  the  integrand  will  be  transformed  into  a 
tractable  form.  From  Bessel  function  theory  (Ref.  13). 


(21) 


e{-[(a2-b2)/4]cos  t}  a  t  (*2  "  ^ 


'a2  -  b2 

(•  l)nIn  - 

\  4 


cos  n£ 


Now  the  function 


(22) 


— .  1  +  2  F  |.|<1 

1  +  a  m*l 
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is  considered.  Substituting 

(23)  z  s  l^eH;  |  n  |  <  1 

in  Eq.  22,  expanding  and  equating  real  terms 


(24) 

Substituting 

(25) 


1  - 


1  +  +  2p  cos  £ 


00 

=  1  +  2  2  (-D’V"  cos  m£ 


m=  1 


a  -  b 


H-  = 


p  |  <  1  ( since  a  and  b  >  0) 


a  +  b 

in  Eq.  24  and  clearing  fractions 
2ab 

(26)  - =1  +  2  V  (-1) 

a2  +  b2  +  (a2  -  b2)  cos  £  m-l 


00 


ml 


la  -  b\m 


la  +  b  i 


cos  m£ 


Substituting  Eq.  21  and  26  for  the  integrand  in  Eq.  20 


(27) 

it  */n 


i2  -  b2\ 


00 


+  2  2  (-DnIn 

n=  1 


i2  -  b2\ 


4  / 


-Icos  n£ 


00 

1  +  2  2  (-U 

m=  1 


a  -  b 


im 


m  i 


_  l  .  e[-(a2+b2)/4] 


i  a  +  b  / 
/a2  -  b2' 


cos  mt, 


d£ 


*  la  -  b\n  /a2  -  b2\l 


Ir 


+  2 


Ir 


n= 1  \ a  +  b 


It  remains  to  express  p(a,  b)  in  terms  of  the  circular  coverage  function. 
Now  the  circular  coverage  function  (Eq.  5)  can  be  expressed  in  the  follow¬ 
ing  forms 


(28) 


p(R,  r)  =  e["(r2/2)1jT  et "  (*2/2)  1  iQ( rt)t  dt 

=  1  -  e['(t2/2^I0(rt)tdt 

-  eHRWl  £  /^\nin(Rr) 

n=l  \  r  / 


1  -  e 


[-(R2+r2)/2] 
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The  latter  two  forms  of  Eq.  28  are  obtained  from  the  previous  forms  by 
repeated  integration  by  parts,  making  use  of  the  relations  (Ref.  13). 

(29)  I.n(t)  =  InW;  <n  an  integer) 


—  [tn  in^)]  *  fnin- 1^) 
dt 


Now  let 


[t-nXn(t>]  -  t-nIn+1(t) 


a  -  b 


a  +  b 


Then  from  Eq.  27,  28,  and  30 
(31)  p(a,  b)  *  1  -  eI_(a2+b2)/4^  Ic 


a2  -  b2 


^  /a-  b 
+  2  £  — r 

n=  1  \  a  +  b 


a2  -  b2 


a  +  b  a  -  b 


b  a  +  b 


From  Eq.  13  and  31 


a  +  b  a  -  b  \  /a-b  a  +  b 
P(a,  b)  *  p  - ,  -  -  p  - ,  — — 

“I  2  2  /  2  2 


a-b  a  +  b 


a+b  a-b 


Therefore  the  probability,  p(a,  b),  of  a  missile  hitting  the  centered 
elliptical  target  with  semiaxes  a  and  b  can  be  computed  from  available 
tables  (Ref.  8,  9,  and  11).  It  will  be  noted  that  if  a  *  b  ■  R,  the  target 
is  a  circle  and  Eq.  31  and  32  reduce  to  Eq.  6,  the  circular  probability 
distribution.  From  Eq.  10  and  13  it  is  obvious  that  the  absolute  value 
of  a  -  b  may  be  used  in  Eq.  31  and  32. 

From  a  geometrical  point  of  view,  the  integral,  p(a,  b),  of  the 
Gaussian  distribution  over  a  centered  ellipse  is  equal  to  the  difference 
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of  the  integrals  of  the  distribution  over  two  associated  offset  circles  as 
given  by  Eq.  32.  A.  centered  ellipse  with  associated  offset  circles  A 
and  B  is  shown  plotted  on  a  circular  probability  grid  in  Fig.  2.  By 
Eq.  32,  the  integral  of  the  distribution  over  the  ellipse  is  equal  to  the 
integral  of  the  distribution  over  the  lune  bounded  by  the  offset  circles. 

A  comparison  of  the  values  of  p(a,  b)  for  the  ellipse  in  Fig.  2,  ob¬ 
tained  by  interpolation  from  the  tables  of  Ref.  8  and  by  count  of  proba¬ 
bility  elements,  is  of  interest.  The  semiaxes  of  the  ellipse  are  a=  1.690 
and  b  =  0.845  standard  deviations,  respectively.  Then  (a  +  b)/Z-  1.268 
and  (a  -  b)/2  *  0.423.  The  values  obtained  by  both  methods  are  given 
in  Table  1.  The  results  show  good  agreement. 


TABLE  1.  Evaluation  of  the  Integral  of  the 
Circular  Gaussian  Distribution  Over  a 
Centered  Ellipse  From  Tables  and  by 
Probability  Element  Count 


Value  of  the  Integral 

Area  of  Integration 

From  Tables  of 
Ref.  8 

By  Probability 
Element  Count 

Offset  circle  A 
p(  1.268,  0.423) 

0.5206 

0.520 

Offset  circle  B 
p(0.423,  1.268) 

0.0397 

0.041 

Lune 

p(  1.268,  0.423)- 
£(0.423,  1.268) 

0.4809 

0.479 

Ellipse 

p(  1.690,  0.845) 

0.4809a 

0.478 

aIntegral  over  ellipse  equal  to  that  over  lune  by 
Eq.  32. 


INTEGRAL  OF  THE  ELLIPTIC  GAUSSIAN  DISTRIBUTION 
OVER  A  CENTERED  CIRCLE 


The  general  form  of  this  integral  is  given  by 


1 

(33)  p(R,  <rx,  fry)  = - 

2  IT  fr^fr  y 

Under  the  transformation 


//  e( ' <It2/2<r*  +y2/2<ry)]dxdy 
x^+y^<R^ 
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FIG.  2.  Centered  Ellipse  With  Associated  Offset  Circles  Plotted  on  a 
Circular  Gaussian  Distribution  Grid. 
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(34) 

x  *  <rxu 

y  *  OyV 

Eq.  33  becomes 

1 

(35)  p(R,  0x>  0y)  *  — 

2tt 

//  e[-(u2+v2)/2]dudv 

(u2/a2  +  v2/b2)<  1 

where 

(36) 

a  =  R/ox 
b  =  R/ay 

Eq.  35  is  similar  to  Eq.  14  and  may  be  evaluated  in  the  same  way.  Hence 


(37)  p(R,  <rx,<ry)  *pl 


*(1+1|.*(1 
2  \(Tx  <Ty  I  2  \<rx 


■R  /  1 


-(1+1 
2\<Tx  <Ty 


INTEGRAL  OF  THE  ELLIPTIC  GAUSSIAN  DISTRIBUTION 
OVER  A  CENTERED  ELLIPSE 


General  forms  of  this  integral  are  given  by 


(38)  p(a,  b,  p,  <rx,  <ry) 

...  1  //  e{-[lA(l-p2)](x2A|-2pxy/^c<ry+y2^)|dxdy 

2ir<rxCy*/l-p2  (x2/a2+y2^>2)<  i 


and 


(39)  p(a,  b,  <rx,  cr  )  =  —  f  f  et_(x2/2,rx+y2/2<ry)1dxdy 

7  2-n  J  J 

where  the  domain  of  integration  is  defined  by 

(40)  Ax2  +  Bxy  +  Cy2  <  D;  B2  -  4 AC  <  0 

These  general  forms  can  always  be  transformed  into  Eq.  14  by  suitable 
linear  transformations.  The  transformations  are  available  in  standard 
texts  in  analytical  geometry  and  are  not  discussed  here. 


INTEGRAL  OF  THE  ELLIPTIC  GAUSSIAN  DISTRIBUTION 
OVER  A  CENTERED  RANDOMLY  ORIENTED  ELLIPSE 

Returning  to  Eq.  33,  p(R,  orx>  o’y)  *s  the  probability  for  given  orx 
and  <Ty  that  a  missile  will  hit  within  a  distance,  R,  of  the  center  of 
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impact.  Hence  p(R,  o-x,  cry)  is  the  probability  distribution  in  R.  If  the 
elliptic  distribution  has  random  orientation,  all  orientations  being 
equally  probable,  the  distribution  with  respect  to  R  remains  the  same, 
but  a  circular  non-Gaussian  distribution  in  hits  is  obtained.  Such  a 
distribution  would  obtain  for  a  randomly  oriented  line  of  delivery  with 
fixed  center  of  impact,  or  for  a  rotating  target  plane  where  the  center 
of  rotation  and  center  of  impact  are  coincident  and  all  orientations  of 
the  target  plane  with  respect  to  any  given  hit  are  equally  probable. 

It  is  of  interest  to  determine  the  probability  density  as  a  function  of 
the  distance,  R,  from  the  center  of  impact.  This  may  be  done  by  ex¬ 
pressing  p(R,  ox.  oy)  in  Eq.  37  as  a  function  of  R  by  means  of  Eq.  5 
and  36,  differentiating  with  respect  to  R  and  making  use  of  Eq.  29. 
After  some  tedious  manipulation 


(41) 

dp .  1  J-d/^dAI+iA^RZjj 

[44  ■  4M 

dR  <rx<Ty 

■  4Wx  °y/ 

is  obtained.  An  alternative  method  of  deriving  Eq.  41  which  will  be 
useful  in  further  developments  is  given  in  the  Appendix. 

If,  as  previously  discussed,  the  elliptic  distribution  has  random 
orientation,  all  orientations  being  equally  probable,  the  probability 
density  with  respect  to  R  as  given  by  Eq.  41  remains  the  same,  but 
becomes  uniformly  distributed  over  the  area  2'irRdR  of  arbitrary  an¬ 
nular  elements  of  radius  R.  Then  the  probability  density  at  the  point 
(R,  0)  or  (x,  y)  is  given  by 


1  dp 

(42) - = 


.[-(l/^Kl/orZ+l/a^R^j 


2irR  dR  2wa-xcTy 


1/1  1  \ 
- r2 

4  \  <r3  <t£  / 


,[-  (1/4)  (lAx  +  V^y)  (x2+y2)  ]  j 


2TT  (Tjj  (T  y 


-(—  -  —  )(*2  +  y2) 

4  \<r2  0-2  / 


Steve  Gaspar,  formerly  of  this  Station,  in  an  internal  memoran¬ 
dum  dated  30  March  1959  also  derived  the  above  probability  density 
function.  His  derivation  differs  from  those  given  in  this  report  in  that 
he  derives  the  probability  density  function  directly  by  considering  the 
probability  density  at  an  arbitrary  point  x,  y  on  the  basis  of  a  centered 
but  randomly  oriented  elliptic  Gaussian  distribution. 
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If  a  centered  elliptic  target  is  allowed  to  have  any  orientation  with  re¬ 
spect  to  an  elliptic  Gaussian  distribution  where  all  orientations  of  the  target 
are  equally  probable,  the  probability  of  a  missile  hitting  the  target  under 
these  assumed  conditions  would  be  the  same  as  if  the  orientation  of  the  tar¬ 
get  were  fixed  and  the  orientation  of  the  distribution  with  respect  to  the  tar- 
get  were  random,  all  orientations  being  equally  probable.  The  probability 
density  of  missile  dispersions  under  this  condition  is  given  by  Eq.  42.  Then 
the  probability  of  a  missile  hitting  the  target  is  given  by 

(43) 


p(a,  b,  crx,  <Ty) 
1 


SI. 


,[-(l/4)(lAx+  lAy  )(x^+y  2)  ]  ! 


2ff(rxo*y  (  X2/a2+y2/b2)<l 


1/  1 


4  \o-|  <r 


-  —  |(x2  +  y2) 


dx  dy 


2"  ft(9)  J-(  l/4)(  l/o&l- l/<r|)t2]j 


ZTTO'xO'y  •'O  •'O 


f 7 

Ja  •'o 


1  j 

1 

1 

.4' 

l<Jx 

4I 

tdtdG 


where  the  second  form  of  Eq.  43  is  obtained  by  means  of  the  transformation 
of  Eq.  15  and  t(0)  is  given  by  Eq.  17.  If  in  Eq.  56  of  the  Appendix,  R  is  re¬ 
placed  by  t(9) 


(44)  £ 


t(0)  /  1  1 

2  \<rx  <rv 


t(0)  / 
2  ,o 


-  £ 


t(  0)  /  l 


2  \c 


2  \ffV 


x  °y 


1  |*t(0)  e[-(l/4)(l/cr^+l/cr^)t2]Ic 

rr  rr_ 


<Tx<Ty  -O 


<rx  <ryl 


tdt 


is  obtained.  Substituting  Eq.  44  in  Eq.  43, 


(45)  p(a,  b,  c rx,  <ry) 


^LL-_LU- 

2  \<rx  Oy/J  2  \<Tx  ®y/  2  \<rx 


due  to  the  periodicity  of  the  integrand. 

Equation  45  gives  the  probability  for  given  <rx  and  <ry  of  a  missile  hit  on 
an  elliptical  target  of  semi-axes  a  and  b  if  center  of  target  and  center  of 
impact  were  coincident  but  the  orientation  of  the  target  with  respect  to  the 
distribution  were  random,  all  orientations  being  equally  probable.  Such  a 
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probability  would  obtain  for  a  missile  hitting  a  rotating  elliptical  target 
where  the  center  of  the  target,  center  of  rotation,  and  center  of  impact 
were  coincident  and  all  orientations  of  the  target  with  respect  to  any 
given  hit  were  equally  probable.  Equation  45  might  find  application, 
for  instance,  in  evaluating  the  probability  of  a  hit  with  a  depth  charge 
on  a  submarine  (approximated  by  an  ellipse)  whose  location  was  known 
but  orientation  unknown.  Computation  of  the  probability  of  a  hit  where 
for  fixed  line  of  delivery  the  lateral  and  range  deviations  are  correlated 
would  involve  a  preliminary  rotational  transformation  of  coordinates  to 
present  the  distribution  in  terms  of  uncorrelated  variables.  Several 
\  rticular  cases  of  Eq.  45  should  be  noted.  If  the  target  is  a  circle, 
then  a  =  b  =  R,  t(0)  *  R  and  Eq.  45  reduces  to  Eq.  37,  the  probability 
that  the  missile  will  hit  within  a  distance  R  of  the  center  of  impact  as 
would  be  expected.  If  <rx  =  (Ty  =  <r,  say,  Eq.  45  becomes 

1  -2ir  /t(  0)  V  /  t(  0)  \ 

(46)  p(a,  b,  <r)s  —  J  jd - ,  ol  -  jdO,  - Id© 

.  e {"[t^6)/2<r^]} de 


e[-(l/2)a2b2/(a2sin20+b2cos20)]d0 


=  1.1  r2Tre|-O/2)(a/<r)2(bA)2/[(a/(r)2sin20+(V<r)2cos20]}d0 

2ir  Jo 

where  substitutions  in  the  integrand  are  made  from  Eq.  6  and  8.  The 
last  form  of  Eq.  46  is  similar  to  that  of  Eq.  18.  Proceeding  in  the 
same  manner,  Eq.  45  reduces  to  Eq.  32,  the  probability  of  the  missile 
hitting  a  centered  elliptical  target  assuming  circular  Gaussian  distri¬ 
bution.  This  is  consistent,  for  with  this  distribution,  orientation  of 
the  target  and/or  distribution  becomes  immaterial.  Finally,  if 
o-x  =  cry  =  <r  and  a  =  b  *  R,  Eq.  45  reduces  to  Eq.  6,  the  probability  of 
the  missile  hitting  a  circular  target  assuming  circular  Gaussian  dis¬ 
tribution. 

Equation  45  may  be  evaluated  by  numerical  integration  using  tabu¬ 
lated  values  of  the  circular  coverage  function  and  suitable  increments 
of  0.  If  Eq.  45  becomes  of  sufficient  importance  to  merit  tabulation, 
it  may  be  non-dimensionalized  as  follows.  Let 

(47)  k  »  b/a 

(48)  i  a  ffy/o-jj* 
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(49) 

Then  Eq.  45  becomes 

(50)  p(a,  b,  <rx,  oy) 


where 


m  *  a/o\ 


2  /*w/2 

Jo  £ 


k(i  +  l)m  k(l  -  l)m 

iKO) - 4-(6) - 


k(i  -  l)m  k(l  +  l)m 

+(0) - +(0) - 


de 


i 


(51) 


m  = 


n/  sin^  0  +  k^  cos^  0 
Equation  51  could  be  tabulated  for  a  range  of  values  of  k,  f,  and  m. 


CONCLUSIONS 

The  integral  of  the  circular  Gaussian  distribution  over  a  centered 
ellipse  may  be  directly  evaluated  from  tables  of  the  circular  coverage 
function.  The  integral  of  the  elliptic  Gaussian  distribution  over  a 
centered  circle  may  be  readily  evaluated  from  tables  of  the  circular 
coverage  function.  The  integral  of  the  elliptic  Gaussian  distribution 
over  a  centered  ellipse  may  be  readily  evaluated  from  tables  of  the 
circular  Gaussian  distribution  if  the  axes  of  the  ellipse  are  parallel  to 
the  axes  of  the  elliptic  equi- probability  contours.  If  the  axes  of  the 
ellipse  are  inclined,  evaluation  can  still  be  made  after  a  preliminary 
transformation  which  may  be  laborious.  Finally  the  probable  value  of 
the  integral  of  the  elliptic  Gaussian  distribution  over  a  centered  but 
randomly  oriented  ellipse  where  all  orientations  are  equally  probable 
may  be  evaluated  by  numerical  integration  of  an  integral  involving  the 
circular  coverage  function. 
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Appendix 

PROBABILITY  DENSITY  OF  A  RANDOMLY  ORIENTED  — 
ELLIPTIC  GAUSSIAN  DISTRIBUTION 

The  probability  density  of  Eq.  38  may  also  be  derived  from  Eq.  37  by  the 
following  method.  We  first  transform  Eq.  37  by  means  of  Eq.  15.  Then 


(52)  p(R,  trx»  0y) 


=  P 


_/i 


2  \o\ 


2  \o\ 


-  P 


2  \(r* 


2  \« 


_ [_  rRe[-( l/4)(l/<r2+l/cr2)t2]  t  rWc[-(l/4)(l/^.lAy)t2cos20]d9dt 

2Trcrx(Ty  •'o 


Substituting 

(53) 

into  Eq.  52 

(54)  p(R,  <rx,  cry) 


^  =  20 


=  P 


2  \(Tv 


2  \<r„ 


-  P 


R/  1 


2  \o\ 


i).R(i+i' 

Jv/  2  \(TX  (Tv 


1  rRe[-(lA>(V^|+lAy)t2]t  /•2ire[-0/4)(l/<r2-l/<r2)t2cosC]dCdt 

V(Tv  *^o  2tt 


4TT<rxO*y 


=  1  fR  e[-(l/4)(lAx+l/«r2)t2]  t  r^J-d/^dAi-lA^cosC]^^ 

•ffO'xO'y  J0  Jo 

due  to  the  periodicity  of  the  integrand.  From  Bessel  function  theory  (Ref.  13) 


(55)  eHl/4)U/«x-  V^y)*2  cost]  m  j 


1/  1 


1 


4W2  ‘  0y/ 
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Substituting  Eq.  55  in  the  integrand  of  Eq.  54  and  integrating 


(56) 


1  fR  e[-(V'4)<l/<r|+V'<ry>«2]ti 


if1  - 1 

4  ^ 


dt 


Differentiating  Eq.  56  with  respect  to  R,  we  obtain  the  probability  den¬ 
sity  as  a  function  of  distance  from  the  center  of  impact.  That  is, 


(57) 


dp 

dR 


1  e[-(l/4Kl/<r2+l/or2)R2]Io 

°x  ^y  ^  ' °j 


R 


From  the  middle  member  of  Eq.  56,  it  follows  that  p(«o,  <rx»  fry)  * 
This  can  also  be  shown  directly  from  the  right  hand  member,  making 
use  of  the  definite  integral  (Ref.  13) 


(58) 


Jf»0O 

r  e(-ap)lc 

o 


(P)dp  = 


•/irrT 


Q  >  1 
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NOMENCLATURE 

a,  b  Axes  of  elliptic  target  area  or  area  of  integration. 

In(t)  Modified  Bessel  function  of  the  first  kind  of  order  n. 

p  Probability  that  a  missile  will  hit  the  target. 

£( R.  r)  Circular  coverage  function.  The  probability  that  a  missile 

will  hit  a  circle  of  radius  R  if  it  is  aimed  at  a  point  a  distance 
r  from  the  center  of  the  circle  and  if  it  is  subject  to  a  circular 
Gaussian  impact-probability  distribution  of  unit  standard 
deviation . 

q(R,  r)  1  -  p(R,  r).  The  probability  that  a  missile  will  not  hit  a  circle 
~  of  radius  R  if  it  is  aimed  at  a  point  a  distance  r  from  the  center 

of  the  circle  and  if  it  is  subject  to  a  circular  Gaussian  impact- 
probability  distribution  of  unit  standard  deviation. 

r  Distance  of  center  of  target  from  center  of  impact. 

R  Radius  of  target. 

t  Radius  vector  coordinate. 

x  Lateral  deviation  of  missile  from  center  of  impact, 
y  Range  deviation  of  missile  from  center  of  impact, 
z  x  +  iy  where  i  =  V-l .  Complex  variable. 

p  Coefficient  of  correlation  between  lateral  and  range  deviations, 
o’  Lateral  and  range  standard  deviations  (used  when  o*x  s  o’y). 
o’x  Lateral  standard  deviation, 
o’y  Range  standard  deviation. 

0  Polar  angle  coordinate. 
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